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ABSTRACT 



We present observations of the optical "moving lines" in spectra of the Galac- 
tic relativistic jet source SS433 spread over a twenty year baseline from 1979 to 
> ■ 1999. The red/blue-shifts of the lines reveal the apparent precession of the jet 

^ ■ axis in SS433, and we present a new determination of the precession parame- 

CN ■ ters based on these data. We investigate the amplitude and nature of time- and 

Q ■ phase-dependent deviations from the kinematic model for the jet precession, in- 

^^ eluding an upper limit on any precessional period derivative ofP<5xlO^^. 

We also dicuss the implications of these results for the origins of the relativistic 



Q_|. jets in SS433. 



Subject headings: binaries - stars: individual (SS433) 



^ ' 1. Introduction 



SS433 is the first known example of a Galactic relativistic jet source, and thus the 
forerunner of modern microquasar astrophysics. The optical spectrum of this object shows 
a number of strong, broad emission lines of the Balmer and Hel series, as well as several 
lines at unusual wavelengths. These latter have been identified as red/blue-shifted Balmer 
and Hel emission from coUimated jets with intrinsic velocities of f ~ 0.26c (Abell & Margon 
1979). Furthermore, the Doppler shifts of these features change with time in a cosinusoidal 
manner, leading to the label of "moving lines" . This behavior is now widely accepted to be a 
symptom of precession of the jet axis in SS433 on a timescale of ~ 164 days (Margon 1984). 

Early studies of the precession in SS433 indicated possible instabilities or drifts in the 
precessional clock (Anderson et al. 1983), which could give considerable insight into the ac- 
cretion processes which must provide the precessional torque. However, Margon & Anderson 
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(1989) reviewed ten years of SS433 timing data and concluded that while significant devia- 
tions from cosinusoidal behavior exist in SS433, the evidence for systematic long-term drifts 
(e.g. precessional period derivative, P) remained inconclusive. 

In this paper, we take the data set considered by Margon & Anderson (1989) and add 
to it more than 50 Doppler shift measurements spread over 10 years, including 9 Doppler 
shifts measured in 1999. Combined, these observations span more than 20 years, and thus 
provide an excellent data set for constraining long-term drifts in SS433's precessional clock. 
We discuss the observations in Section 2. In Section 3, we present analyses of the entire 
data set in the context of the "kinematic model" for SS433's precessing jets. In Section 4, 
we discuss the results of these analyses, and in Section 5 we present our conclusions. 



2. Observations 

The primary observations used here are optical spectroscopic observations of SS433, from 
a wide range of telescopes and instruments (see Margon & Anderson (1989) and references 
therein for details). The net result of these observations, spread over the period from June 
1978 to July 1992, is the measurement of 433 Doppler shifts for the "receding jet" {zi) and 
482 Doppler shifts for the "approaching jet" (^2) for the optical moving lines in SS433. 

We obtained further spectra of SS433 in July, 1999 using the Hartung-Boothroyd Ob- 
servatory (HBO) 24-inch telescope and optical CCD spectrograph. We used a 600 lines/mm 
grating and 6" slit providing a resolution of i? ~ 800 (6 A/pix). We present a typical 
spectrum in Figure 1. 

We determined the Doppler shift of each HBO spectrum using only the moving Ha lines, 
and we did so by fitting a Gaussian profile to the red and blue components separately. Note 
that the profiles of the moving lines are broad, time-variable, and often asymmetric. This is 
due to the time overlap of multiple discrete emission components, commonly referred to as 
"bullets" (Vermeulen et al. 1993), with typical lifetimes of ~ 3 days. These systematic devia- 
tions introduce a relatively large uncertainty in the Doppler shift determination. Based upon 
examination of many spectra, we find that the typical full-width at half-maximum (FWHM) 
is Az ~ 0.003, and we adopt this as our uncertainty in the Doppler shift determination o"z. 
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3. Analysis 

3.1. The Kinematic Model 

Throughout our analysis of these data, we adopt the "kinematic model" for the moving 
lines, which assumes that the changing Doppler shifts arise from the precession of the jet 
axis in SS433. The simplest form of the kinematic model takes into account five components: 
the jet velocity P = -; the jet angle from the precessional axis 9, the inclination angle of 
the system with respect to the observer's line of sight i, the precession period P, and the 
epoch of zero precessional phase to- The period and zero-phase epoch combine to give the 
precessional phase = ^^. The resulting Doppler shifts obey the equation 

Zi^2 = 1 — 7[1 ± P sin^ sini cos0 ± j3 cos 6 cosi] 

where 7 = (1 — P^)^^'"^. SS433 exhibits "nodding" of the jets on a ~ 6.5-day period (Katz et 
al. 1982) due to the ~ 13-day binary motion of SS433 (Crampton et al. 1980), which are not 
accounted for in this model. However, the effects of this nodding are essentially negligible for 
long timescale studies of the jets such as ours. We further mitigate the impact of nodding by 
applying a 7-day boxcar smoothing filter to the individual Doppler shifts determined above. 
We then used chi-squared minimization to find the best-fit parameters for the kinematic 
model (Table 1). The resulting model fit is plotted versus time along with the data and 
residuals in Figures 2-3 for zi and Z2. We plot the same model fit, data, and residuals versus 
precessional phase in Figure 4. 

The resulting fit has a chi-squared residual per degree of freedom of xt = 8-9) indicating 
the presence of statistically significant residuals. However, we can still use this fit to estimate 
uncertainties in the kinematic model parameters as follows. First, we scale all of the a^ 
values by v^8^, so that the residuals have xl-ftx — 1-0; essentially by fiat. We then take 
the uncertainties to be the range of a model parameter which introduces a total change of 
Ax'fix = 1-0. We also report these values in Table 1. This rescaling approach for deriving 
the model parameter uncertainties is statistically valid in a strict sense only if the residuals 
are consistent with Gaussian noise and are not correlated with any model parameters in a 
systamtic way. If so, then the residuals would simply indicate that we have ignored one 
or more sources of noise in the system when estimating the uncertainties in the individual 
Doppler shifts. As we show below, this is largely true, though we see some evidence of 
small (but statistically-significant) systematic deviations from the kinematic model. Thus, 
the uncertainties in the model parameters in Table 1 are likely to be good, but not perfect, 
statistical estimates. For the remainder of the paper, we adopt the best-fit model parameters 
presented in Table 1. 
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3.2. Doppler Shift Residuals 

One obvious feature of Figures 2-4 is tliat tlie residuals to the model fit greatly exceed 
the uncertainties in the Doppler shift determinations (as also shown by the large value of xt 
above). We also notice no obvious trend in the residuals versus time as would be expected 
for systematic timing effects, such as a constant precessional period time derivative P. Such 
large, apparently random residuals have been noticed in previous timing studies of SS433 
(Anderson et al. (1983); Margon & Anderson (1989)). 



3.2.1. Correlations in Residuals 

Previous studies have also noticed that the velocity residuals in SS433 show a pattern 
of correlation between zi and Z2 (Margon & Anderson 1989). Specifically, when we plot the 
residuals of Zi{ohs) — Zi{mod) versus Z2{obs) — Z2{mod), we find that most of the points lie 
in the second and fourth quadrants (Figure 5). In other words, when the absolute value of 
Zi is greater than expected, the absolute vlaue of Z2 is also greater than expected, and vice 
versa. The number of data points with zi and Z2 residuals in quadrants 2 and 4 is 271 ± 16, 
while in quadrants 1 and 3 the number is 110 ± 10 - a > 8o" difference. 

The linear correlation coefficient between the residuals is r = —0.69 ±0.02. We estimate 
the uncertainty from a Monte Carlo simulation as follows. We take the 381 pairs of zi and 
Z2 residuals and add to each a random number drawn from a Gaussian distribution with 
mean of zero and a standard deviation of 0.003 - the typical uncertainty in the Doppler 
shift measurements. We then calculate the correlation coefficient of the resulting simulated 
distribution. We repeat this procedure 1000 times and then take the standard deviation in 
the correlation coefficient as the uncertainty above, ar = ±0.02. 

This correlation pattern could have several physical sources. The effect considered most 
commonly in previous studies (e.g. Margon & Anderson (1989)) is that of phase noise in 
the precessional motion, with strict symmetry between zi and Z2. As the jet precessional 
phase either lags or leads the model ephemeris, the projected velocity amplitudes of the 
jets on the observers line of sight will either exceed or fall short of the model prediction. 
Another possible physical explanation is modulation of the velocity amplitude - "/9-noise" - 
in a system which otherwise follows the 5-parameter kinematic model ideally (e.g. Milgrom 
et al. (1982)). 
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3.2.2. Phase-Dependence of Residuals 

Another factor which could impact the Zi/z2 residual correlation are phase-dependent 
residuals to the kinematic model. Margon & Anderson (1989) found no evidence for such 
phase-dependence in their data. In analyzing this data set, we divided the data into 10 
evenly-spaced phase intervals and calculated the average and standard deviation of the resid- 
ual Doppler shifts in each bin (Table 2). None of the average residuals from the kinematic 
model is as large as the rms deviation for residuals in that phase bin, indicating that phase- 
dependent deviations from the kinematic model do not dominate the residuals. Furthermore, 
if we compare the average residuals of the individual phase bins to the rms scatter of all the 
phase bins, none of them are more than 2a outliers. Thus, the amplitude of any average 
deviation from the model velocity does not seem to be phase-dependent. 

However, when we calculate the uncertainty in the average residual for each phase bin, 
equal to the standard deviation in the residuals divided by the square-root of the number 
of points, we find that the deviations are in fact statistically significant (Table 2). That is, 
while the scatter around the kinematic model in any given phase bin is not dominated by 
systematic deviations from the model, such deviations are present in the data set. The nature 
of these systematic deviations are not clearly determined. Figures 2-3 show that the velocity 
residuals show some correlation on timescales of weeks or months. Thus, sparse sampling of 
the velocities combined with such correlations in the residuals could be one explanation for 
the apparent systematic deviations. 

We used the average residuals from the kinematic model as correction factors for the 
quadrant analysis presented above and in Figure 5, taking the observations and subtracting 
both the kinematic model and the average deviation for all points within A0 = ±0.05 cycles 
of each data point. As a result, the number of ,21,^2 residual pairs in quadrants 2 and 4 
does decrease, but only to 261 (with 120 pairs in quadrants 1 and 3). Thus the correlation 
between zi and Z2 residuals remains highly statistically significant even after correction for 
the systematic deviations. 

We also note that the relative phase-independence of the residuals raises questions re- 
garding the nature of the residuals. As can be seen from the equation for Doppler shifts in 
the kinematic model, every parameter which could be time-variable {l3,P,i,6,) either feeds 
directly into the phase </> or is multiplied by cos 0. Thus, noise in these terms or in itself 
should result in a cosinusoidal modulation in the RMS of the Doppler shift residuals with 0. 
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3.2.3. The Phase Noise Model 

As mentioned above, the most commonly-invoked physical model for the velocity resid- 
uals in SS433 is "phase jitter" in the jet precession. Since the precession phase affects both 
jets similarly, it naturally explains the correlation between the Zi and Z2 residuals. If such 
jitter can occur over timescales of weeks or months, it can also explain the long-term residual 
correlations evident in Figures 2-3. 

We analyzed this SS433 data set following the example of Margon & Anderson (1989), 
determining phase errors from the velocity residuals above. We simply defined the phase 
error to be the phase difference between the actual phase of the observation given its epoch 
and the kinematic model parameters in Table 2 and the closest model point with the same 
observed velocity. As can be seen in Figure 4, some observed velocity amplitudes exceed 
the maximum model velocity amplitude, and such points were dropped from this analysis. 
We then divided the data set into 10-day intervals and calculated the average and standard 
deviation of the phase errors from all phase measurements doing that interval (including 
both Zi and ^2)- For 10-day intervals with only 1 phase measurement we have no estimate 
of the standard deviation, and thus dropped such intervals from the analysis. We plot the 
resulting phase noise measurements in Figure 6. We repeated this same analysis using a 
30-day interval for averaging, with the results shown in Figure 7. 

We note that while there are occasional trends in the residuals on timescales of several 
hundred days, no obvious trend is apparent over the full time span in either panel of Figure 
7. The 1999 data are marginally inconsistent with zero phase residual (at the 2. 80" level for 
one of the two data points in Figure 6b). However, it is clear that this phase residual is less 
than many prior apprently secular deviations from the kinematic model in Figures 6-7. If 
we assume that some period derivative is present in SS433 over the span of our observations, 
these secular deviations could mask its effects up to A0 ~ 0.05 cycles. Given the span of our 
observations, this corresponds to an upper limit on the period derivative of P < 5 x 10~^. 



3.2.4. The Velocity Amplitude Noise Model 

As mentioned above, an alternate physical explanation for the velocity residuals in 
Figures 2-4 is noise in the intrinsic velocity of the jets. To investigate this possibility further, 
we calculated the intrinsic jet velocity necessary to match each observed Doppler shift, given 
6, i, to, and p from the best-fit parameter set in Table 1. We plot the corresponding values 
for (3 = - versus time in Figure 8 and versus precessional phase in Figure 9. 

The average value of 13 we find is 0.254 with a standard deviation of 0.024. Given 
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507 independent measurements of (3 in this way, we arrive at an average value of (3ave = 
0.254 ± 0.0011. While this value is only ~ 4% lower than the value given in Table 1, the 
difference is statistically significant at the 7.9a level. This may indicate that "noise" in the 
Doppler shifts may in fact be impacting the parameter estimates for the kinematic model in 
a systematic way, as discussed in Section 3.1. 



4. Discussion 

4.1. Phase Noise 

As noted above, the upper limit on precessional period derivative of P < 5 x 10~^ 
shows that there is no large long-term drift in the precessional timing properties of SS433. 
The presence of jitter in the system implies some "torque noise" in the process driving the 
precession, according to the phase noise model. However, if this were the case, that noise 
must average out over timescales of ~ 20 years. We can also see from Figure 7 that there 
are fairly large phase deviations of A0 ~ 0.1 cycles over timescales as short as ~ 10 days. 
This implies that the torque noise At has a maximum relative amplitude of at least 

Ar^ax (A0/At) 



1-^ 



T {d(j)ldt) 

~ 1.6 

Thus, the variation in torque can in fact exceed the time-averaged torque driving the preces- 
sion. This may be a problem for certain physical models of the precession and timing noise 
in SS433. 

Finally, we note that the phase noise model is incapable of producing the observed 
Doppler shifts which exceed the maximum amplitude predicted by the kinematic model. 
We have considered the possibility that the phase noise itself causes the x-squared fitting 
procedure used to determine the model parameters to systematically underestimate the true 
jet velocity, and thus "undershoot" the maxima. However, Monte Carlo simulations of data 
sets with higher true velocities and phase noise identical to that observed here fail to produce 
such undershooting. Therefore, we conclude that phase noise model cannot reproduce the 
observed Doppler shift residuals near the maximum projected velocities. 



4.2. Velocity Noise 

The alternate "/9-noise" model, on the other hand, can clearly explain the excess velocity 
at the extrema (and any other precessional phase) by the changing jet velocity amplitude. 
Such a model also has a physical basis, given recent advances in the modeling of relativistic jet 
production. Meier et al. (2001) dicuss a scenario where such jets are launched by a magnetic 
accretion disk instability around a black hole (or other compact object). Variations in the 
accretion flow onto the compact object (i.e. M, intrinsic magnetic field, etc.) can alter the 
radius at which the magnetic field saturates and the jet is launched, and thus the jet velocity. 
The relation between jet velocity and launch radius for a non-rotating black hole follows: 



where f3 = -, and Rg is the gravitational radius of the black hole (one-half of the Schwarzschild 
radius) . 



4.3. Jitter models and phase-dependence of residuals 

It is also interesting to view these model in light of the apparent lack of phase-dependence 
in the Doppler shift residuals noted above. By differentiating the equation for Doppler 
shifts in the kinematic model (eqn. 1) with respect to the potentially time- varying model 
components {f3, i, 6, and (p) we can see the relative phase-dependence of Doppler shift 
residuals on deviations in each term. In the phase noise model, we would expect the following 
dependence: 

Az = {'-) (3 sin 9 sin i sin 0) A0 
^0.1 sin0 A0 

Thus, we would expect the amplitude of the Doppler shift residuals to be sinusoidally mod- 
ulated with phase, in apparent contradiction with our analyses above. 

For the "/3-noise" model, we have: 

Az = (7 sin^ sini cos0 — 7 cos^ cosi) A(3 
- (0.35 cos - 0.2) A/3 
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(with the approximation that d^jdfi ~ 0). Again, we have a modulation of the Doppler 
shift residual amplitude dominated by a term varying cosinusoidally with respect to phase. 

For variations in the angle between the jet axis and the precessional axis, 6*, we have: 



I\z = (7 /? cos^ sini cos0 + 7 /5 sin^ cosi) IS.Q 
~ (0.24 cos + 0.09) M 

again dominated by a cos term. 

Finally, for variations in the system inclination angle, i, we have: 

A2; = (7 /3 sin^ cosi cos0 + 7/9 cos^ sini) Ai 
~ (0.02 cos + 0.24) A? 

Interestingly, the Doppler shift residual amplitudes for "i-noise" would be dominated by a 
constant term, with only a small dependence on phase. Thus, this is the only parameter in 
the kinematic model for which variations causing the Doppler shift residuals are consistent 
with their observed phase-independence. Unfortunately, we know of no physical model for 
such variability at this time. Furthermore, based on the mass estimates for the compact 
object and companion star (~ IOMq total), the known 13.5-day binary period, and the 
presence of deviations as large as 0.01 rad/day in inclination angle, the change in rotational 
energy would require average powers of -E ~ 10'^ Lem^ which seems implausible. 

One other possible explanation is that the Doppler shift residuals are due to variations in 
both precessional phase and one (or more) of the other model parameters. In that case, the 
residual amplitude would depend on a sum of cos and sin terms which could potentially 
smooth out any phase-dependence in the residuals. 



5. Conclusions 

We have presented observations of the Doppler-shifted optical moving lines in SS433 
spanning over 20 years. We draw the following conclusions based on the data: 

• We find parameters for the kinematic model for the jet precession which are similar to 
those found by previous authors (e.g. Margon & Anderson (1989)). 
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We find a strong correlation between residuals to the models fits for the two jets Zi 
and Z2, with a linear correlation coefficient of r = —0.69 ± 0.02. 

We find that the residuals to the kinematic model fit are not dominated by systematic 
phase-dependent deviations from the model. However, systematic phase-dependent 
deviations from the kinematic model are seen in the data set at a low level. 

If we adopt a "phase noise" model for the velocity residuals, we find correlated devia- 
tions over timescales of months to years, but no long-term trend over the full data set. 
We place a limit on the precessional period derivative of P < 5 x 10^^. 

Noise in any single parameter of the kinematic model seems unable to explain the 
observed phase-independence of the velocity residuals in SS433. However, variations 
in both phase and one of the other parameters would vary as the weighted sum of cos 
and sin terms, which could smooth out any phase- dependence of the Doppler shift 
residuals. 
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Table 1. Best-fit parameters for tlie kinematic jet precession model 

Parameter f5 6 i P to 

(deg) (deg) (days) (TJD) 

Value 0.2647 20.92 78.05 162.375 3563.23 

Uncertainty ±0.0008 ±0.08 ±0.05 ±0.011 ±0.11 



Table 2. Average residuals in Doppler shift versus phase 



Phase Interval 


< Azi > 


A^i 


Azr 


<AZ2 > 


AZ2 


AZ2 


(cycles) 




RMS 


Uncert. 




RMS 


Uncert. 


0.0-0.1 


-0.0035 


0.0074 


0.0010 


0.0040 


0.0067 


0.0008 


0.1-0.2 


0.0017 


0.0082 


0.0011 


-0.0024 


0.0083 


0.0010 


0.2-0.3 


-0.0005 


0.0101 


0.0014 


-0.0016 


0.0093 


0.0014 


0.3-0.4 


0.0035 


0.0076 


0.0011 


0.0003 


0.0028 


0.0004 


0.4-0.5 


-0.0011 


0.0130 


0.0032 


-0.0008 


0.0067 


0.0013 


0.5-0.6 


-0.0052 


0.0116 


0.0016 


0.0031 


0.0083 


0.0010 


0.6-0.7 


0.0039 


0.0152 


0.0028 


-0.0008 


0.0135 


0.0025 


0.7-0.8 


0.0015 


0.0092 


0.0013 


-0.0001 


0.0077 


0.0010 


0.8-0.9 


0.0026 


0.0104 


0.0013 


-0.0018 


0.0079 


0.0010 


0.9-1.0 


-0.0035 


0.0088 


0.0013 


0.0007 


0.0063 


0.0008 
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Fig. 1. — Typical spectrum of SS433 taken from the Hartung-Boothroyd Observatory 24-inch 
telescope. Note the clear appearance of the moving Balmer lines from the jets. 
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Fig. 2. — Doppler shift data, model fit, and residuals versus time based on the best-fit param- 
eters in Table 1 for zi, the "receding jet" . The vertical extent of the plotting symbols shows 
the typical ilo" uncertainty in the Doppler shift. 
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Fig. 3. — Doppler shift data, model fit, and residuals versus time based on the best-fit pa- 
rameters in Table 1 for Z2, the "approaching jet" . The vertical extent of the plotting symbols 
shows the typical ±lo" uncertainty in the Doppler shift. 
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Fig. 4. — Doppler shift data, model fit, and residuals versus precessional phase based on the 
best-fit parameters in Table 1 for: (top) zi, the "receding jet" ; (bottom) Z2, the "approaching" 
jet. The vertical extent of the plotting symbols shows the typical ±lcr uncertainty in the 
Doppler shifts. 
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Fig. 5. — Doppler shift residuals from the kinematic model. Note the clear anti- correlation 
between Z\ and Z2 residuals. 
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Fig. 6. — Phase residuals deduced from velocity residuals to the kinematic model, averaged 
over 10-day intervals. 
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Fig. 7. — Phase residuals deduced from velocity residuals to the kinematic model, averaged 
over 30-day intervals. 
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Fig. 8. — Jet velocity amplitude required to match the observed Doppler shifts as a function 
of time. 
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Fig. 9. — Jet velocity amplitude required to match the observed Doppler shifts as a function 
of precessional phase. 



